On the Density Dependent Nuclear Matter Compressibility 
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In the present work we apply a quantum hadrodynamic effective model in the mean-field approx- 
imation to the description of neutron stars. We consider an adjustable derivative-coupling model 
and study the parameter influence on the dynamics of the system by analyzing the full range of 
values they can take. We establish a set of parameters which define a specific model that is able 
to describe phenomenological properties such as the effective nucleon mass at saturation as well as 
global static properties of neutron stars (mass and radius). If one uses observational data to fix 
the maximum mass for neutron stars by a specific model, we are able to predict the compression 
modulus of nuclear matter K — 257, 2MeV. 
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INTRODUCTION 



Effective theories that implement hadronic dynamics 
are considered the most efficient way to describe infinite 
nuclear matter because they allow to work in a specific 
energy scale while ignoring other degrees of freedom of 
the system. In the present work nuclear matter is consid- 
ered at low temperature and low chemical potential, thus 
the relevant degrees of freedom are hadrons. Typically 
two pathways are explored to get insight in hadronic mat- 
ter properties, on the one side at high temperatures and 
low densities, heavy ion collisions, on the other hand ex- 
ploring neutron stars at low temperatures and extreme 
densities With the increase of temperature in the 
first case, or density in the second, the quarks deconfine 
and chiral symmetry is effectively restored. Although 
this feature can be relevant for the study of neutron star 
matter, the existing models able to describe it [^SU] are 
very complex and do not allow a straightforward analysis 
of more complex functional dependence of the couplings, 
that is exactly the intention of this work. In this article 
we present a generalized model theory for hadron dynam- 
ics in a nuclear matter environment, supposed to exist in 
neutron stars, and discuss the density dependent nuclear 
matter compressibility and its response to the parame- 
ters of the theory. In an initial step the model theory 
is calibrated for densities where from the experimental 
point of view phenomenology is well under control. In a 
second step, the undetermined parameters are partially 
fixed by a symmetry scheme or varied freely. The princi- 
pal feature of our model theory is that it unites proper- 
ties of a variety of approaches known in the literature ([B[ 
and references therein) at the cost of additional parame- 
ters. Nevertheless, as the following discussion will show, 
only a straight range for those parameters is physically 



'Electronic address: dexheimer@th.physik.uni-frankfurt.de 



meaningful, so that the increase of degrees of freedom is 
considerably reduced by phenomenology. 

Further we make profit from the fact that small varia- 
tions in the parameters create considerable differences in 
the results. For example a change in the binding energy 
of nuclear matter of 0.25MeV implies a decrease by 0.5 
solar masses for the maximum mass as predicted by a 
model for a specific family of neutron stars as shown in 
more detail below. For that reason a precise knowledge 
of input parameters is crucial in order to make predic- 
tions relevant for observation. In the present work the 
input parameters are: the binding energy and the asym- 
metry coefficient of nuclear matter, the effective mass of 
the baryons and the compression modulus at the nuclear 
matter saturation point. The compression modulus may 
be understood as a measure for the difficulty to remove 
the system from equilibrium or for the affinity to return 
to equilibrium. 

Until now, compressibility is one of the parameters de- 
termined with major uncertainty. One of the goals of 
the present work is to show the influence of the onset of 
hyperon presence and how this might be useful in pin- 
ning down the uncertainty in the compression modulus 
using observational evidences. To this end we calculate 
the compressibility for a large range of densities including 
the usually determined value at saturation density. 



II. THE MODEL 

The basis for our calculations is a Quantum Hadro- 
dynamics motivated model including beside the nucleons 
also the hyperons which define the baryon octet (p, n, 
A, E+, E°, £ _ , S° and H — ). The nuclear medium effect 
is taken into account by a non-linear coupling involving 
the sigma field with adjustable parameters A, and 7, 
which tune the intensity of the meson baryon couplings 
(for details see 0, The considered mesons are the 
usual scalar a and the isoscalar vector meson u> but also 
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the isovector vector meson p. Additionally the lepton 
degrees of freedom (e and //) establish or complement 
charge neutrality, depending on the density. The La- 
grangian density reads: 

B 

-{M B - gl B <r)]ipB + ^ tpiiij^ - m e )ip e 
i 

+ \(d^a - mla 2 ) - + ^m*w„w" 

-\p^ v + \™ 2 p p^, (1) 

where the index B denotes a specific baryonic species and 
the index I a specific leptonic species. Here the adjustable 
coupling strengths are related with the usual coupling 
constants through 

9*oB = X*\b9<jB, 
s£b = X*fm9uB, 

9* p b = X* iB 9 P b, (2) 

and the adjustable coefficients xIb that drive the dynam- 
ical influence of the hadronic medium are 

for i = A, (3 or 7 (see eq. [2j). These parameters are 
constrained in order to reproduce phenomenology, in the 
present case an effective mass at saturation between 0.7 
and 0.8 nucleon masses and a compression modulus be- 
tween 200MeV" and SOOMeV. The coupling constants 
9a n, 9ujN and g P N of the nucleons are chosen in a way to 
have zero pressure and a determined energy at satura- 
tion (calculated assuming a binding energy of —16MeV 
and a density of 0.17/m -3 at saturation and a bare nu- 
cleon mass of 939MeV) and to reproduce an asymmetry 
coefficient of 32.5MeV. Since no stringent experimental 
confirmation with respect to meson hyperon couplings 
is available and also possible nuclear medium effects on 
that couplings are an open question, one way to reduce 
these ambiguities is to use the SU (6) flavor spin symme- 
try scheme as in ref. as a guide. 

.9(tS = 2/3g CTjv 5<tH = 1/3.9^ g a A = 2/3^^ 
5wS = 2/3g W Ar OluS = l/3g^N g^A = 2/ig^N (4) 

9p£, = 2g P N g P ~ = g P N g P \ = 0. 

As already indicated before, one advantage of our 
model is that it is general in the sense, that certain 
choices of parameters A, /3, 7, permit to recover the princi- 
pal models cited in the literature. For A = (3 = 7 = the 
Walecka model [|fis reproduced, for A = 1 and f3 = 7 = 
the ZM1 model [3] is recovered, for A = /3 — 7 = 1 the 
ZM3 model j§] results and A = (3 — 7 — > 00 describes the 
exponential model [10]. Beside the possibility to create 



models using intermediate values of the parameters, we 
can group the parameters according to the combinations 
in which they are being varied: 

1. A varies with j3 = 7 = 0, 

2. A varies with A = (3 = 7, 

3. A, /3 and 7 vary independently. 

For the first coupling scheme we find acceptable values 
for the effective mass and the compression modulus at 
saturation for 0.06 < A < 0.19. For the coupling scheme 
2 there are no acceptable values, which shows that the 
sigma field influence on the vector mesons is not "univer- 
sal". For the last coupling scheme there are acceptable 
values for 0.00 < A < 20.22 and 0.00 < (3 < 1.35, where 
both parameters correlate [1, 0| . A variety of combina- 
tions have shown us that only a considerable restricted 
choice of parameter combinations is possible in order to 
reproduce phenomenological meaningful quantities. The 
7 parameter cannot be related to symmetric nuclear mat- 
ter saturation properties as it relates to the p meson and 
consequently to the isospin asymmetry. 

Depending on the density, the constituent particles of 
the system change. The hyperons decay into other hy- 
perons so that baryon number conservation and chemical 
equilibrium have to be taken into account. Note, that 
the chemical potential is the criterion which indicates 
the threshold where heavier baryons start to contribute 
to the chemical composition of nuclear matter. Further 
charge neutrality shall hold, according to findings in the 
literature, that a spherically symmetric charged neutron 
star is not stable [ll|. To solve the system analytically 
and in a self-consistent way the mean-field approxima- 
tion is used, since for higher densities the fluctuations of 
the meson fields are negligible (see for instance ref. [§]). 

III. COMPRESSIBILITY FUNCTION 

The compressibility function of the star as a function 
of density is determined from its thermodynamical defi- 
nition (proportional to the derivative of pressure with re- 
spect to volume) 0, [f| and written as the first derivative 
of pressure P (with respect to the total baryon density), 
the second derivative of the energy density e, or using the 
definition of baryon chemical potential fi n = de/dp: 

where /i„ is equal to the chemical potential for neutrons. 
The chemical potential for different baryon species can 
be calculated from: 

fi B = yJk 2 FB + (M*) 2 + (g^V^o + (g* pB ) P03l3B ■ (6) 

For each density the chemical equilibrium equations de- 
termine the chemical potential for each species which de- 
termines its Fermi energy. Once the particle's chemical 
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potential is less than the Fermi energy the system fa- 
vors the presence of the respective particle. Here the 
chemical potential of each species is written as a func- 
tion of the effective coupling strengths, and isospin is 
explicitly included. Isospin effects and charge conserva- 
tion together with the nontrivial response of the effective 
mass to the presence of the hadronic environment man- 
ifest themselves in the phenomenon that hyperons com- 
posing the star matter as a function of density do not 
appear simply according to their vacuum mass hierar- 
chy. Fig. Q] (a) shows this behavior in the compressibility 
function as a function of baryon density. 

In agreement with relativistic formulations this model 
respects the causal limit ( i.e., the speed of sound in the 
medium is smaller than the speed of light for any den- 
sity) , so from the relation of compressibility and baryon 
chemical potential the speed of sound in the medium is 



G) 3 



K 



and the following inequality is established 
K < 9fi n . 



(7) 



(8) 



l°g,n e <- (S /cm ) 



In Fig. Q] (a) this fact is indicated by the limiting dotted 
line in the respective plot of the compressibility function. 



IV. NEUTRON STARS 



FIG. 1: Compressibility as a function of baryon density (a) 
and mass of the star as a function of central energy density 
(b) for A = (Walecka model). 



FIG. 2: Observed masses of known pulsars in units of solar 
masses. 



The equation of state obtained from the La- 
grangian density (U is then plugged into the Tolman- 
Oppenheimer-Volkoff equation, which adds gravitational 
effects to nuclear matter properties. For each of the pa- 
rameter sets, a family of neutron stars is obtained and 
compared with observational data. Then it is possible to 
analyze the behavior of the 7 parameter in the neutron 
star mass and radius so that the best parameter set may 
be chosen. 

As can be seen in the Fig. El within the observa- 
tional error bars the mass of almost all example pulsars 
have a value approximately 1.4M Q . However, lately pul- 
sars with higher masses (the largest one with a mass of 
2.1 ± O.2Af ) have been observed. Since there is no cri- 
terion as to which model (QHD, quark matter or others) 
is the most adequate starting point in order to derive the 
equation of state, we do not consider the largest possible 
mass as a relevant reference for the present model and 
focus on finding the model that agrees with nuclear phe- 
nomenology and predict the largest possible neutron star 
mass (~ 1.5Mq). The radius of the model stars range 
between 10 km and 14 km, which because of uncertain- 
ties in obtaining stringent values from observation, does 
not permit to put further constraints on the parameters. 

The model with all the parameters equal to zero, which 
corresponds to the Walecka model extended by the inclu- 
sion of hyperons from the baryon octet and the meson p 
predicts a larger maximum mass 1.9M© (Fig. 21(b)) but 
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FIG. 3: Compressibility as a function of baryon density (a) 
and mass of the star as a function of central energy density 
(b) for A — > oo. 



FIG. 4: Compressibility as a function of baryon density (a) 
and mass of the star as a function of central energy density 
(b) for A, P -> oo. 



suffers a drawback because of a considerable small ef- 
fective mass M* ~ 0.5M and an amplified compression 
modulus K ~ 500MeU. 

An increase in A, while keeping j3 e 7 equal to zero, de- 
creases the compressibility function (Fig. 0(a)) because 
the star can only support a lower mass against gravity 
(Fig. [3](b)). This behavior happens for values A < 1.75. 
After that, the system saturates, meaning that it does 
no respond to further changes in the parameter. The 
range with acceptable nuclear phenomenology is between 
0.06 < A < 0.19 and corresponds to maximum mass from 
1.46 to 1.22 solar masses. 

Upon increasing A and (3 together, keeping 7 equal to 
zero, the compressibility function (Fig. [4] (a)) and the 
neutron star mass (Fig. |4] (b)) decrease. A variety of 
combinations for A and (3 (not shown in this article, but 
elsewhere 0) have shown us that (3 and A have a similar 
effect on the macroscopic properties, but an increase in 
j3 shifts the onset of saturation towards larger values for 
A, which saturates for A > 16 and j3 > 10. The range 
with acceptable nuclear phenomenology ranges between 



0.00 < A < 20.22 and 0.00 < /3 < 1.35. 

As the third option a variation of 7, which tunes the 
isospin contribution in nuclear matter, also (slightly) 
modifies the compressibility function (Fig. [5] (a)) and 
neutron star mass (Fig. [5](b)), respectively. Depending 
on the isospin sign of the particle (for example positive 
for proton and negative for neutron) , the contribution of 
the term related to 7 can be positive or negative making 
the net effect small in comparison with A and contri- 
butions. 



V. CONCLUSION 

In the present article we calculated the density depen- 
dent compressibility function for a generalized model for 
nuclear matter. Usually, compressibility is determined 
only at saturation density, however, in order to charac- 
terize states in the final phase of a neutron star evolu- 
tion, where quasi equilibrium conditions are encountered, 
a large density range for the compressibility gives details 
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FIG. 5: Compressibility as a function of baryon density (a) 
and mass of the star as a function of central energy density 
(b) for A, j3, 7 — > oo (exponential model). 



FIG. 6: Compressibility as a function of baryon density (a) 
and mass of the star as a function of central energy density 
(b) for A = 0.06, /3 = 0, 7 = -1. 



relevant for possible processes. 

With increasing density, hyperons come into play and 
soften the equation of state. The couplings have shown 
to be of crucial influence in the conditions which deter- 
mine the sequence in which hyperons appear. Because 
the present model mimics medium effects by the depen- 
dence of the couplings on the sigma field, this nonlin- 
earity is responsible for the relevant hyperon hierarchy, 
which does not reflect the sequence of threshold by sim- 
ply considering the masses. 

Among the variety of possible models, the one in agree- 
ment with nuclear phenomenology and astrophysical ob- 
servation is defined by the parameter set A = 0.06, 
(3 = 0.00 and 7 = —1.00. This model has an effec- 
tive nucleon mass of 0.7Mn, a compression modulus of 
257.2MeV (for comparison see refs. [H, O, El) and 
a maximum neutron star mass of 1.48M Q as shown in 
Fig. [6l This choice allow us to calibrate the model that 
can now be used, in a second step, to calculate many 
other neutron star properties such as radius, redshift, 
maximum rotational frequency, moment of inertia, cool- 



ing properties, among others. 

At first sight the value of the maximum neutron star 
mass does not include the largest observed masses. But 
analyzing Fig. [2] one may reason that within the error 
bars the masses of most of the pulsars are compatible 
with ~ 1.4 solar masses. Stars with larger masses may 
well contain a new kind of matter, such as quark matter, 
quark hadron mixtures or other exotic phases [13 [13 . 

With the present work we established a model which 
is in agreement with nuclear phenomenology and obser- 
vational data. Using the maximum mass as a selection 
criterion, determined by varying the model parameters 
A, /3 and 7, we were able to predict a compressibility of 
K = 257, 2MeV from astrophysical arguments, in fairly 
good agreement with the limits given by other analysis. 
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